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Abstract —In a centralized radio access network (RAN), the 
signals from multiple radio access points (RAPs) are processed 
centrally in a data center. Centralized RAN enables advanced 
interference coordination strategies while leveraging the elastic 
provisioning of data processing resources. It is particularly well 
suited for dense deployments, such as within a large building 
where the RAPs are connected via fibre and many cells are 
underutilized. This paper considers the computational require¬ 
ments of centralized RAN with the goal of illuminating the 
benefits of pooling computational resources. A new analytical 
framework is proposed for quantifying the computational load 
associated with the centralized processing of uplink signals in 
the presence of block Rayleigh fading, distance-dependent path- 
loss, and fractional power control. Several new performance 
metrics are defined, including computational outage probability, 
outage complexity, computational gain, computational diversity, 
and the complexity-rate tradeoff. The validity of the analytical 
framework is confirmed by comparing it numerically with a 
simulator compliant with the 3GPP LTE standard. Using the 
developed metrics, it is shown that centralizing the computing 
resources provides a higher net throughput per computational 
resource as compared to local processing. 

Index Terms —Computational complexity, computational out¬ 
age, computational diversity, turbo-decoding, mobile networks, 
3GPP LTE 

I. Introduction 

Compared with today’s networks, future mobile networks 
will need to cope with a dramatic increase in the density 
and demand of users. Several novel technologies have been 
proposed to meet the needs of future systems, including 
massive MIMO, millimeter wave signaling, very densely 
deployed small-cell networks, and centralized radio access 
networks (RANs) m, El. The centralized RAN concept is 
illustrated in Fig. [T] A centralized RAN system consists of 
multiple radio access points (RAPs), which perform only 
simple processing of the input signal, including sampling 
and optionally performing an FFT operation, a fronthaul or 
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backhaul connection, a transport network, and a central data 
center. The data center manages a virtual base-station (BS) 
pool composed of a scalable number of virtual machines, each 
representing one or more BSs and performing the majority of 
the baseband processing. 

To date, research on centralized RAN has focused on the 
RAN functional split, the applicability of joint processing, 
system performance, and fronthaul requirements. Recently, 
the implementation on cloud-computing platforms attracted 
more interest El: cloud-computing platforms allow resources 
to be virtualized, improve resource elasticity, and allow mobile 
networks to be implemented in software. However, the appli¬ 
cation of centralized RAN to cloud-computing platforms (also 
referred to as Cloud-RAN) requires a better understanding 
of the interplay between available computing resources and 
the required communications performance. As we will show, 
the amount of computational resources and the degree to 
which they are centralized have a direct impact on how the 
communication performance is effected by the communication 
channel properties and resource allocation policies. 

A. Data processing complexity in mobile networks 

Each mobile network is characterized by a well-dehned 
set of functionality, timing requirements, and protocols. This 
imposes very precise requirements on the operation of each 
RAP, including data processing requirements in order to 
maintain real-time properties. For instance, 3GPP LTE dehnes 
a set of modulation and coding schemes (MCSs), each with 
a prescribed number of information bits per transport block 
(TB). Since the choice of MCSs depends on the signal-to- 
noise ratio (SNR) and the computational requirements de¬ 
pend on the number of received information bits, there is 
an inherent relationship between the channel quality and the 
computational requirements. Eurthermore, 3GPP LTE specifies 
strict processing deadlines that are necessary to guarantee 
given latency constraints and prevent unnecessary hybrid-ARQ 
retransmissions. Due to randomness in the channel and user 
demand, the data processing requirements are typically found 
by resorting to extensive simulation studies. 

Because the local processing capability of each small cell 
is anticipated to be far less than that of a macro cell, 
the dimensioning of data processing resources needs to be 
revisited. Eurthermore, very dense networks are subject to 
dramatic temporal and spatial traffic fluctuations such that 
many small cells may be strongly underutilized. It is therefore 
not economically viable to equip a small cell based on 
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Fig. 1. Typical centralized RAN architecture 

peak data processing requirements, yet under-dimensioning the 
computational resources limits its capabilities. This dilemma 
motivates the concept of centralized RAN, which enables 
computational load balancing across multiple BSs to avoid 
peak-provisioning. The goal is to dimension the computational 
assets of the network to exploit the computational diversity 
present in a pool of resources, while still guaranteeing the 
communication service requirements of each individual cell. 
The utilization of processing resources should be maximized 
for a given service guarantee. These goals require a framework 
that enables the prediction of the required data processing 
resources under a given performance constraint. 

B. Related work 

The major benefits of centralized RAN have been unveiled 
in m. Key among these benefits is computational load bal¬ 
ancing, which exploits the variable number of users processed 
within the data center (a consequence of the spatial and 
temporal traffic fluctuations) and variations of the processing 
requirements per user (due to the channel variability). In ID, 
Bhaumik et al. use the the OpenAir 0 LTE implementa¬ 
tion and real-world measurements to predict that centralized 
RAN can save up to 22 % of the data processing resources. 
Furthermore, the paper considers statistical guarantees that 
sufficient computational resources are provided with very high 
probability. However, these results do not characterize the 
source of fluctuations, they do not allow the quantification 
of the centralized RAN gains depending on the variation of 
channel parameters, and they do not allow for active shaping 
of the data processing load. In 0, the impact of the decoder 
complexity on the data processing requirements is evaluated 
numerically. Furthermore, the impact of computational con¬ 
straints is investigated on the overall system performance, 
and it is shown that taking these constraints into account 
by the scheduler allows the system to recover much of the 
performance loss caused by imposing a computational budget. 
However, the work was largely based on simulation, and does 
not provide a formal analytical framework for the evaluation 
of centralized RAN and its computational requirements. 

Implementation options for centralized RAN have been 
investigated in Q, which considers flexible centralization over 
heterogeneous backhaul with possibly high latency. In IS), the 
signal processing requirements of this system are investigated. 
Among others, different centralization options are analyzed 


with respect to their required backhaul capacity, latency, and 
challenges for the signal processing. A preferred functional 
split is to centralize the forward error correction coding and 
decoding (possibly on commodity hardware), while all other 
lower layer functionalities are executed at the RAP. In 0, 
a virtual BS pool based on commodity hardware and Linux 
operating system is investigated. It is shown that such a system 
can provide the required real time and throughput guarantees. 

Evaluating the complexity requirements of a centralized 
RAN system requires an analysis of the corresponding base¬ 
band processing. However, the evaluation of computational 
complexity in mobile networks has been primarily focused on 
the power consumption of receivers. So far, mainly energy- 
limited devices have been investigated where complexity and 
power consumption are of particular interest, e.g. in ||9l where 
ultra-wideband receivers have been analyzed. The receiver 
complexity is affected by different components such as the 
decoding and detection process as well as Altering process. 
Hong and Stark estimated in ifTOl the power consumption and 
decoding performance of prominent decoders for codes such 
as turbo, Hadamard, block, and convolutional codes. In in, 
a thorough analysis of the sphere decoding and detection has 
been performed with particular focus on complexity reduction 
in multiple antenna systems. Similarly, ma analyzed low- 
density parity check (LDPC) codes and possible complexity 
reductions. The Alter complexity of the receiver has been 
investigated in in the context of software-defined radio. 
In software-defined radio, it is important to reduce complexity 
particularly at the lower layers. At higher layers, scheduling 
incurs a major part of the complexity, which has been formally 
analyzed in m, which derives expressions quantifying the 
complexity required to schedule a successful transmission. In 
Elias and Jalden investigate in the trade-off between rate, 
reliability (expressed by outage probability), and complexity 
in the high-SNR regime. 

Complexity is equally important from an infrastructure point 
of view. In na, the network capacity is investigated using 
game theory and distributed algorithms, where the complex¬ 
ity is distributed across multiple nodes. In El, Howard et 
al. provide a formal model for quantifying power consumption 
and derive the critical distance at which the transmit power 
consumption equals the decoder power consumption. This 
critical distance has been quantified for particular decoder 
implementations. In lITSl . Grover et al. derive a comprehensive 
framework for analytically modeling the power consumption 
of a decoder and introduce an optimization problem whose 
main objective is to minimize the system-wide power con¬ 
sumption including both transmission and decoding power. 
The optimal rate allocation to minimize the system-wide 
power consumption is strictly lower than Shannon capacity. In 
this paper, we leverage the model proposed in lIT^ to evaluate 
the data processing requirements of a centralized RAN system. 

C. Contribution and Outline 

In this paper, we introduce a framework for modeling, an¬ 
alyzing the performance, and determining the data processing 
requirements of centralized RAN systems, which applies for 
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both macro and small cells. The model is used to quantify the 
computational load associated with the centralized processing 
of uplink signals in the presence of block Rayleigh fading, 
distance-dependent path-loss, and fractional power control. 
Several new performance metrics are defined, including com¬ 
putational outage probability, outage complexity, computa¬ 
tional gain, computational diversity, and the complexity-rate 
tradeoff. The validity of the analytical framework is confirmed 
by comparing it numerically with a simulator compliant with 
the 3GPP LTE standard. Using the developed metrics, the 
advantages of centralized RAN over conventional (distributed) 
RAN are illuminated. 

The remainder of this paper is organized as follows. Section 
iniintroduces the system model, presents the complexity model, 
and defines the new performance metrics. In Section [Till an 
analytical framework for computational complexity is intro¬ 
duced, culminating in tractable expressions of the proposed 
complexity metrics suitable for cellular uplinks characterized 
by block Rayleigh fading and fractional power control. Section 
IIVI performs a numerical verification of the proposed frame¬ 
work. Section [V] presents and discusses results obtained with 
our framework. Finally, the paper concludes in Section |VT] 

II. Complexity model and Pereormance metrics 

Consider a RAP group containing RAPs, where each 
RAP corresponds to a single BS and whose signals are jointly 
processed in a virtual BS pool. Let 7 indicate the instantaneous 
SNR for a particular link and 7 = E[ 7 ] indicate the average 
SNR. In the following, it is assumed that the channel gain 
remains fixed for the duration of one TB, but varies indepen¬ 
dently from TB to TB, which corresponds to a block-fading 
model. For each TB, the channel is conditionally subject to 
additive white Gaussian noise (AWGN). When the SNR of a 
channel is 7 , its capacity is log 2 (l + 7 ) per channel use lfT9l , 
which provides an upper bound on the rate. In practice, the 
rate is selected for a given 7 from a set of MCSs in such a way 
that a constraint on outage probability is satisfied. Because no 
finite-length code satisfies the capacity limit with equality, the 
rate is selected to be below the corresponding capacity limit 
according to some margin. 

A. Link Adaptation 

In order to perform link adaptation, the transmitter can 
choose from Nr different MCSs, with Rk denoting the rate 
of the MCS, k G {1,..., W/j}. As in LTE, we assume that 
each MCS corresponds to setting the rate of a turbo code and 
selecting from among three modulations: QPSK, 16-QAM, 
and 64-QAM. See ll20ll . 11211 for details of how the turbo¬ 
code rate and modulation format is defined for each MCS. 
Per the standard, the k*^ MCS segments each TB into Ck 
code blocks (CBs), each of which conveys Dk information 
bits. As we advocate later, the selection of the MCS should 
take into account the complexity of the turbo decoder, and in 
particular, the number of expected iterations required to meet 
the outage constraint is the key factor in making the selection. 

Suppose that the turbo decoder is run with a certain number 
of maximum iterations Lmax- Let 7 ^ indicate the minimum 


SNR for which the k*^ MCS satisfies the outage constraint 
after the Tmax-th iteration, on average. The value of 7 ^ for 
each MCS can be obtained as follows. Simulations are used 
to obtain code block error rate (CBLER) curves for each 
possible MCS and for an arbitrarily large number of maximum 
iterations (here, set to Lmax = 8 ). For the MCS, 
is selected to be the value of SNR for which the transport 
block error rate (TBLER) satisfies a particular constraint for 
the channel outage ^channel (a transport block error occurs when 
any of the code blocks in a transport block fails). 

It is well known that the complexity of a turbo decoder 
depends on how much the channel SNR exceeds the capacity 
limit. If a turbo code of rate Rk is transmitted over a channel 
whose capacity is sufficiently higher than Rk, then typically 
only a few iterations are required. Thus, to manage complexity, 
an SNR margin A 7 may be used. With such a margin, the SNR 
must be at least A 7 above in order to select the MCS. 
The rate selection for channel SNR 7 and SNR margin A 7 
can be represented by the function 

(Ri if^<7f 

r(7,A7) = (1) 

[Rnu if7wH<^- 


B. Complexity Model 

The computational effort required to decode a turbo code 
is linear in the number of iterations and in the number of 
(information) bits. Thus a reasonable metric for computational 
effort is the bit-iteration. For a given SNR 7 and SNR margin 
A 7 , let the random variable Lrij, A 7 ) be the number of iter¬ 
ations required to decode a particular CB, and E A 7 )] 

be the average number of iterations required to decode a CB. 
The expected decoding complexity (averaged over the number 
of iterations required per CB) expressed in bit-iterations per 
channel use (pcu) can be evaluated as 


^ , DkCkE[Lrh,Ay)] 

C( 7 , A 7 ) =--- 


( 2 ) 


where S^e is the number of channel uses or resource elements 
(REs) required to convey the TB. The expected number 
of decoding iterations can be found using the same set of 
simulation curves used to determine the MCS thresholds 7 ^. 
For a given rate r ( 7 , A 7 ), the simulation results will show 
the CBLER for each iteration up to the Lmax-th iteration. 
These CBLER curves can be interpreted as the probability 
mass function (pmf) of the number of iterations at 7 , and 
from the pmf coefficients, the average number of iterations 
is computed. 

As an example, we have computed the complexity using (|2]) 
for an example system. For this and all of the examples in this 
paper, the choice of MCS assumes the use of the normal cyclic 
prefix and a 10 MHz bandwidth corresponding to 50 resource 
blocks (RBs). However, since up to 5 RBs must be reserved 
for the physical uplink control channel and the number of 
RBs must be a multiple of 2, 3, or 5 Il20l , 11211 , the transmitted 
signal occupies 45 RBs. Since there are 12 information-bearing 
SC-FDMA symbols per sub-frame and 12 subcarriers per RB, 
it follows that S're = 45 X 12 X 12 = 6480 . The MCS selection 
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is implemented according to ([TJ with the threshold found for 
^max = 8 iterations, an outage constraint ^channel = 0 . 1 , and 
no excess margin (A 7 = 0). This empirically determined 
complexity is shown in Fig. |2] by the blue curve, while the 
black curve shows the complexity predicted by the complexity 
model that will be discussed next. 

To enable a deeper analysis of complexity issues, it is 
essential to build an accurate yet wieldy complexity model. 
The groundwork for such a model is provided by m, which 
accurately predicts the power consumption of a decoder for 
a given code rate, coding scheme, channel, and decoder type. 
According to ESI Eq. (4) and (9)], the expected decoding 
complexity under a given instantaneous SNR 7 and SNR 
margin A 7 can be modeled as 


C(7,A7) = 


?’(7: A 7 ) 

l0g2 (C - 1) 

r h, A 7 ) 

log2 (C - 1) 

r il, A 7 ) 

log2 (C - 1) . 

-21og2 [I (7,A7)] 


l0g2 

(^channel) C ^ 

. Pi^,Ay)K' C 

logio (^channel) C ^ 

l0g2 

[ P{j,Ay)K' C 

l0g2 

( ) 

\ -^(^channeOC J 


+ 


(3) 


where C is a parameter of the model related to the connectivity 
of the decoder when represented as a graph. 


I ( 7 , A 7 ) 


AT ( f channel) 


log 2 (l + 7 ) - r( 7 , A 7 ), 
K' 

^^g 10 ( ^channel) 


(4) 

(5) 


K' is a parameter of the model, and Cchannei is a constraint 
on the channel outage probability. For the sake of simplicity 
Ar(echannei) IS denoted by K in the remainder of the text. 

When evaluating 0, the function r( 7 , A 7 ) is defined by 
EJ. The relationship between each threshold 7 ^ and the 
corresponding rate Rk in 0 can be modeled as 

1 + ( 6 ) 



where j/ is a parameter that models the gap between the 
capacity at 7 ^ and the SNR for the actual code to meet the 
performance objective at rate Rk- The values of 
could be found for each MCS by statistically fitting the 
empirically observed complexity to the complexity predicted 
by the model. However, it is more convenient to have values 
for these parameters that are common to all the MCSs. 

The values of (C, K', v) that provide the best fit to the em¬ 
pirical complexity (blue curve) of Fig. |2] are ( = 6, K' = 0.2, 
and ly = 0.2 dB. Using these values with 0 provides the 
predicted complexity shown by the black curve in Fig. |2] As 
can be seen, these two curves are quite similar, demonstrating 
the effectiveness of the model. We note that a more accurate 
model could be achieved by determining a different set of 
(C, K', v) for each MCS, but this significantly complicates 
the analysis, as will be subsequently evident. Fig. further 
shows a highly volatile behavior, which results from the link- 
adaptation process. When the system operates at the same 
MCS (rate), a decoder requires fewer iterations as the SNR 



Fig. 2. Complexity as a function of the SNR obtained both through 
simulations (blue) and based on (black curve), when the maximum number 
of iterations used to decode a CB is limited to eight. 


increases. However when link adaptation is performed, the 
system requires more complexity each time it moves towards 
an higher MCS, since it operates closer to capacity, which is 
the cause of the peaky behavior in Fig. |2] 


C. Performance Metrics 

In this subsection, several performance metrics are defined, 
which are used throughout the paper to evaluate the perfor¬ 
mance of a centralized RAN system and highlight its benefits 
over a distributed implementation. 

7) Computational outage probability: Fig. |2] shows that the 
required computational complexity to decode a TB may differ 
significantly depending on the SNR. Furthermore, when the 
system is distributed, each BS will be provisioned with a 
limited amount of computational resources Cmax- By contrast, 
in a centralized system the virtual BS pool will be equipped 
with a total limited amount of computational resources, which 
scales linearly with the number of BSs. If the total required 
computational effort exceeds the computational resources be¬ 
fore a decoding deadline is reached, a computational outage 
will occur. From the operator’s perspective, a computational 
outage is no different than a channel outage, because in both 
cases the information is not made available to the end user. 

More formally, the computational outage probability is 
defined by 


f^comp ( Ac, Umax ) 


■ N, 

^C( 7 j,A:) > AcCtnax 


(7) 


where 7 ^ denotes the instantaneous SNR for the 7-th RAP. 
When a computational outage occurs, the centralized processor 
is unable to finish the decoding of at least one uplink user. 
However, when this condition occurs it is still possible that 
some of the uplink transport blocks can be decoded, even 
though at least one of the blocks has failed. 

2) Outage complexity: The outage complexity Cout(e)Ac) 
is the minimum amount of computational resources required 
to satisfy an outage constraint on the per-cell computational 
outage probability fcomp (Ac). More formally, the outage com¬ 
plexity is defined as 

Ac) — arg min i 6comp(Ac, Umax) t^comp (Ac) r. ( 8 ) 


The constraint of the computational outage probability of the 
system ecomp(Ac) is normalized to a per-cell basis. Assuming 
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independent processing with a non-priority-based scheduler, 
the constraint for a centralized system is related to the con¬ 
straint for a distributed system with value fcomp ( 1 ) = e by 

ecomp(iVe) = 1 - (1 - . (9) 


3) Computational gain and diversity: Using the previous 
definition of computational outage, the centralization gain of 
a centralized RAN mobile network is determined. In particular, 
the impact on the performance by centralizing communication¬ 
processing resources can be measured in two ways: 1 ) by con¬ 
sidering the reduction of the required computational resources; 
2 ) by examining the improved outage behavior. 

The computational gain quantihes the reduction of com¬ 
putational resources made possible by the centralization of 
processing. When RAPs are jointly processed, the compu¬ 
tational gain is 


9comp (-^c) 


A^cCout(e, 1) 
Comic N^) 


( 10 ) 


The asymptotic computation gain is found by letting the 
number of RAPs go to infinity. 


5comp “ 5comp (A^c) • (U) 

^ A'c^OO 


An alternative metric is the computational diversity, which 
quantifies the rate of computational improvement as more 
RAPs are jointly processed. This is found by determining the 
rate at which the computational outage probability is improved 
as a function of Nc, i. e., the computational diversity is 


t:^comp i^c'} — 


^logio ^comp (ATcCmax) 

ON, 


( 12 ) 


This is the downward slope of the curve that relates compu¬ 
tational outage probability to the number of jointly processed 
RAPs. The larger this negative slope, the faster the maximum 
computational gain (fTTl i is achieved. In a mobile network, the 
computational diversity provides an indication for the required 
number of centralized BSs per data-center. 

We note that the slope is at its maximum magnitude when 
Nc = 1, and thus for the remainder of the paper when we 
discuss computational diversity, we evaluate the derivative 
of (fT2l l at Nc = 1. Under this condition, we can drop the 
dependence on Nc and simply express dcomp {Nc) by dcomp- 

4) Average achievable rate: Given a fixed average SNR 7 , 
and a hxed SNR margin A 7 , the rate associated to a given 
user depends upon the MCS selected based on the channel 
quality. Due to the random effect of fading, the code rate 
selected per user is aleatory. One of the conventional ways to 
measure system performance involves evaluating the average 
code rate. Given the distribution of the instantaneous SNR 7 , 
the expected code rate can be computed as follows 


Nr 

{R (A 7 )} = ^ P { 7 « < 7 < 7 «+,} r ( 7 f, A 7 ) 




Nr 

= E r {it A7) ( 13 ) 


k=l 


where F^iy) is the cumulative density function (CDF) of the 
instantaneous SNR 7 . 


5) Complexity-Rate Tradeoff: Focusing only on outage 
complexity ignores the influence of rate. An alternative metric 
should consider the relationship between the rate and the 
SNR margin A 7 . As the margin A 7 increases, the complexity 
decreases, but so does the rate since the system operates farther 
away from capacity. Therefore, the complexity-rate tradeoff 
(CRT) can be quantified, and it measures how much additional 
complexity is required in order to further improve the rate. In 
particular, the complexity-rate tradeoff can be evaluated as 


fcomp(Ac; A 7 ) —lim 


{i?(A7')} 


A-y' 


dAff 


dComicNc) 


dAff 


-1 

(14) 


III. A Framework for Complexity Analysis 


In this section, we provide an analytical framework for 
computing the previously introduced performance metrics. We 
begin in Section IIII-AI by assuming that the user has a fixed 
average SNR 7 , which is true when the uplink is fully power 
controlled. We continue in Section lTlI-BI bv considering the use 
of partial power control with a user that is randomly positioned 
in the cell. For ease of exposition, we start with Nc = 1 and 
extend the results in Section UlI-CI to the case of arbitrary Nc- 


A. Local Processing with Rayleigh Fading 

Consider a single user with fixed average SNR 7 , and that 
this user is processed locally (Nc = 1). The instantaneous 
value of 7 is a random variable, which, for the sake of 
tractability, we assume is an exponential random variable 
corresponding to the Rayleigh fading modefl The probability 
density function (PDF) of such a 7 is given by: 

A ( 7 ) = =exp ,7 > 0. (15) 

1) Expected complexity: 

Theorem 1: The expected value of the decoding complexity 
conditioned on the average SNR 7 can be evaluated by 


Nr 


k=l 


log2 (C - 1 ) 


[/i (7^+1, 7 ^) + h (7^+1, 7 ^)] ( 16 ) 


where the functions Ip (•, •) and I 2 (•, •) are given respectively 
by (1^ and (l43T l. 

Proof: The derivation is given in Appendix lAl ■ 

2 ) Variance of the complexity: 

Theorem 2: The variance of the computational complexity 
conditioned on the average SNR 7 can be evaluated by 


L 


n 2 


logs (C - 1) 

Ci {it-l’Tk) ~ ^-y\l{Ci} 

where E.y|^{Ci} is evaluated using (fThl) and R (•, •) is 

by “ 


Proof: The derivation is given in the Appendix iBl 


(17) 

given 


^We note that other distributions for 7 may be used, more accurately 
modeling the effects of interference, shadowing, and signal propagation, but 
this comes at the cost of reduced tractability. Most such distributions will 
require a simulation or numerical integration to compute the performance 
metrics. 
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3) Outage complexity: In the case of an isolated user pro¬ 
cessed locally, the outage complexity in a distributed system 
is the value Cout that satisfies the following expression 

P{C. > Cout(e, 1)} < e (18) 


which can be evaluated once the CDF of the complexity is 
known. 

The CDF of the complexity can be written as follows 

Pc(C,hr) = P{C. <C,hr} (19) 


where Cthi is a complexity threshold chosen based on the qual¬ 
ity of the channel and the computational resources available. 
Using the notion that the complexity monotonically decreases 
in 7 in the interval (fT^ can be replaced by 

^ Nr 

Pc(Cthr) = P {7fc.min < 7 < 7^1 } (20) 

P|7>7l|^ 


where 

Ik.mm = max 


7f,min 


( 21 ) 


AR = 2 “^*r'°S 2 (C l)/[2r(7.A7)] /C - ^22) 

V KC 

Finally, (l20l i can be written as follows 

, Nr 

Pc(Cthr) = 1 _ p ( R\ [-^7 ( 7 ^ 1 ) ~ P 7 (7fc,min)] (23) 
^ ^7 (7l ) 


where F^{-) is the CDF of the instantaneous SNR with PDF 
given by (flST l. 

4) Example: Fig. [3 shows a comparison between analyt¬ 
ical (solid lines) and simulation (dots) results, in terms of: 
a) expected complexity; b) variance of the complexity; c) 
outage complexity when the constraint on the computational 
outage is set to be e = 0 . 1 ; and d) expected achievable-rate. 
The simulations are obtained as follows. During each trial, 
the instantaneous SNR yi is determined for each uplink by 
drawing an exponentially distributed random variable. For the 
given SNR, the rate is determined according to the function 
r{'ji,Ay). Next, the complexity C{yi,Ay) is determined 
using ®. Once a sufficient number of trials is run, the average 
and variance of C ( 7 ^, A 7 ) are found along with the average 
achievable rate, while Cout(e, 1 ) is evaluated such that (l 8 ]l is 
satisfied for = 1 , given a specific constraint on both the 
computational and channel outage. For the following results, 
one million trials are run for each of Nr MCS schemes and 
for several representative values of SNR margin A 7 . In doing 
this experiment, the value of Nn is varied in an effort to gain 
some insight into the required number of MCSs schemes. 

All four sub-figures show both the dependency on the 
number of MCS schemes Nr and the SNR margin A 7 . 
The number of MCS schemes Nr determines how well the 
capacity curve is sampled, while the SNR margin A 7 provides 
a direct indication of how close the code and decoder operate 
to Shannon capacity. This SNR margin decreases when more 
turbo-decoder iterations are allowed. The value of 7 ^ for each 
MCS is equally-spaced with the values of its adjacent MCSs. 


The first and the last MCS have assigned values, respectively 
7 f = -6.4 dB and 7 ^^ = 17.6 dB (based on 3GPP LTE). 
The average SNR is fixed to 7 = 10 dB. Fig. [3 shows a 
good agreement between the analysis and the simulations for 
all values of A 7 and Nn. However, as the number of MCS 
schemes decreases, the mismatch between the analysis and 
simulations increases for the variance of the complexity and 
this is imputable to the non-linearity of log(l + x) and the 
linearization method adopted to derive (fTTl l. 


B. Effect of Path Loss and Power Control 

Now consider the case that the user is subject to a fractional 
power control policy and is located randomly within the cell. 
For the sake of tractability, we assume a circular cell of unit 
radius and that the user is placed according to a uniform 
distribution. A more complicated model for the placement of 
the user may be used, such as one that accounts for the actual 
shape of the Voronoi region. While more accurate, this comes 
at the cost of tractability and will typically require simulation 
or numerical integration. Let Tj indicate the BS and its 
location. Each BS Yi serves only one user, which is located 
at Xi at distance \Yi — Xi\ from the BS. Under a fractional 
power control policy, Xi will transmit using power 

= Po\Y, - X,\^^ (24) 

where Pq is a reference power (typically taken to be the power 
received at unit distance from the transmitter), rj is the path- 
loss exponent, and s ,0 < s < 1 , is the compensation factor 
for fractional power control. The average SNR received by Y) 
for the mobile is given by 

(25) 


where 7 ^^} is the average reference SNR, measured assuming 
unit-distance transmission. 

In this case, it is not possible to determine a closed-form 
expression or a suitable approximation for both the expected 
and the variance of the complexity. However, knowing that the 
PDE of r, = I A, - y,| is 


2a; , 0 < a; < 1 

0 , otherwise 


(26) 


the expected complexity E{Ci} and variance Var{Ci} can be 
obtained by performing a numerical integration over r of 
(fThl l for the expected complexity and (fTTl l for the complexity 
variance. 

Theorem 3: Under this scenario, the CDE of the instanta¬ 
neous SNR can be computed in closed form and it yields 


J^7(7) = 1- 


(i) 


-2/Ml-s)] 


77(1 - s) 


7(1 - s) ’ 7ud 


)-r 


7(1 - s)’ 


(27) 


where Fj-, •) is the incomplete gamma function ll 22 i . 

Proof: The derivation is given in the Appendix fC] ■ 
The outage complexity can be determined in this case 
analytically, since the CDF of the decoding complexity can 
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Number of MCS schemes Nr 
(b) Variance of the complexity 




Number of MCS schemes Nr Number of MCS schemes Nr 

(c) Outage complexity when e = 0.1 (d) Expected achievable rate 

Fig. 3. Expected complexity, variance of the complexity, outage capacity, and expected achievable rate as function of the number of MCS schemes used, 
for three different values of A 7 and block Rayleigh fading. The average SNR is 7 = 10 dB. The constraint on the computational outage is e = 0.1. Solid 
lines are obtained analytically, while the dots are obtained through simulations. 


be evaluated by using dZTl i into (|23]) . Even if dTTT i cannot be 
expressed in closed form in this case, well-known algorithms 
such as bisection, Newton, or gradient-search algorithm can 
be applied. Furthermore, the average achievable rate can be 
determined by substituting dZTl) into (fOl) . 


C. Centralized Processing 

Using the previously derived expressions for the expectation 
and the variance of the complexity, in this section the analyt¬ 
ical framework is expanded in order to be able to compute 
through an approximation the maximum decoding complexity 
when iVc BSs are centralized. More specifically, the central 
limit theorem (Lindeberg/Levy) 12^ is used, which yields 


Cout(e,iVe) = (2(1 - - 1) +E{C,} 

y ^ V 

inv. norm. CDF 

(28) 

where is the inverse error function, which has the 

following series representation ll 22 l 


k=0 


where vo=l and 


2fc+l 


Vk 


k-1 

771=0 


'^m'^k—l—m 

{m + l)(2m -I- 1) 


(29) 


(30) 


The two sub-hgures of Fig. |4] show respectively the absolute 
and the relative outage complexity which ensure the per-cell 


outage constraint e = 0.1. The normalization of the relative 
outage complexity is done over the decentralized case. In 
both sub-figures the dependency to the number of centralized 
BSs and the SNR margin A 7 is shown. As expected, the 
normalized outage complexity decreases as the number of 
centralized BSs or the SNR margin increases. Fig. 0] shows 
curves obtained both analytically (solid lines) and through 
simulations (dots) using the 3GPP FTE parametrization intro¬ 
duced in Section III-BI The simulated results are obtained as 
follows. At each trial, the instantaneous SNR 7 ^ is determined 
by drawing a random variable exponentially distributed. Once 
the correct MCS is found, the complexity Cnorm is determined 
using ©. Once a sufficient number of trials is run, Cout(e, Nc) 
is evaluated such that (|7|i is satisfied, given a specific constraint 
on both the computational and channel outage. For the follow¬ 
ing results, hundred-thousand trials are run for each number 
of centralized BSs and for each SNR margin A 7 used. The 
path-loss exponent is set to rj — 2, the average SNR at unit 
distance is 7 uji = 0 dB and the compensation factor is fixed to 
s = 0.1. The notches on the right side of each sub-figure show 
the behavior as N^. 00 . Both sub-hgures show a good match 

between analysis and simulations. Furthermore, the right sub- 
hgure emphasizes the beneht of using a centralized RAN 
rather than a distributed system. 

IV. Numerical verification 
In order to evaluate the validity of the analytical framework 
derived in Sec. mil in this section numerical results are 
provided based on the 3GPP FTE standard and a 3GPP FTE 
decoder. The numerical results are compared against those 
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(a) Absolute outage complexity 



(a) Rayleigh fading 



(b) Relative outage complexity 


Fig. 4. Outage complexity to ensure the per-cell outage constraint e = 0.1. 
Solid lines are obtained through analysis, while the dots ai‘e the results of 
simulations. The notches on the right side of each sub-figure show the behavior 
as Nc oo. 



(b) Distance dependent path-loss model 
Fig. 5. Outage complexity to ensure per-cell outage constraint e = 0.1 as 
function of the number of RAPs, whose signals are centrally processed. Solid 
lines are evaluated analytically, while dots are obtained through simulations 
using one million trials. The notches on the right side of each sub-figure show 
the behavior as Nc —yoo. 


obtained through the previously derived framework. LTE uses 
adaptive modulation and coding based on turbo codes with 
overall 27 distinct MCSs (TV/j = 27). Each MCS is identified 
by an MCS index, I^cs = {0, ...,26}, and is characterized 
by a different combination of code rate and modulation 
format Eol, ED. Three kinds of modulation are used: QPSK 
(0 < /mcs < 10), 16-QAM (11 < /mcs < 20), and 64-QAM 
(21 < /mcs < 26). When a TB is larger than 6144 bits, it 
is segmented into multiple CBs. Each CB is separately turbo 
encoded, and all Ck CBs in the TB must be correctly decoded 
for the TB to be correct. Let Ccb be the probability that a CB 
is in a channel outage. It follows that Cchannei = 1 — (1 — 
is the probability that the TB is in a channel outage. 

Eor a given value of SNR 7 ^, the MCS is selected to satisfy 
echannei < echannei where ^channel = 0 . 1 , which is a typical value 
for an LTE network. When complexity is a concern, an SNR 
margin A 7 can be applied, which corresponds to the case 
when the MCS is selected such that the outage constraint is 
met after a specific number of decoder iterations. 

Assume all the available resources per cell are allocated to 
a single users. Eurthermore, a full-buffer scenario is assumed: 
each of the RAPs, whose signals are processed at the 
data center, always serves a user over the complete band¬ 
width. However, 3GPP LTE allows for multi-user opportunistic 
scheduling such that one TB may carry information from more 
than one user. When resources are shared among multiple 


users more diversity and therefore a lower outage complexity 
per cell can be achieved compared to the case when all 
the resources are monopolized by a single user. Therefore, 
this assumption represents a worst-case scenario regarding 
computational complexity, since it does not allow for any 
computational diversity within a given cell. 

A block Rayleigh fading channel is considered, although 
in practice a system experiences more channel diversity due 
to multi-path propagation, and the computational complexity 
varies from user to user and from CB to CB. This scenario, 
in an extreme case, leads the users to experience a sufficient 
channel diversity for a close-to-AWGN channel, which results 
in a lower computational gain and diversity. 

The outage complexity Cout(e,Ac) is obtained through a 
simulation campaign. During each trial, the SNR 7 ^ is drawn 
from an exponential distribution when we do not consider 
the path loss, or from a distribution whose CDE is given by 
(l27l i when the path loss is taken into consideration. The MCS 
scheme for a given 7 ^ is determined according to the value 
of A 7 chosen. Each of the Ck CBs in the TB is marked as 
being in a channel outage with probability Ccb, which can be 
precomputed. If any of the CB are in outage, then the entire 
TB is declared to be in an outage. If the CB is in an outage, 
then Lr ( 7 , A 7 ) = L^ax, which is the maximum number of 
attempted iterations. Otherwise, Lr ( 7 , A 7 ) is determined by 
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drawing a random variable distributed according to the pdf of 
Lj. ( 7 , A 7 ), which can be precomputed by tracking the error- 
rate as a function of the number of iterations. Once the decod¬ 
ing complexity is computed, finally the outage complexity is 
evaluated by (l 8 ]l. The outage complexity is obtained by using 
one million trials when evaluated through simulations. 

Fig.|5]shows the outage complexity to ensure per-cell outage 
constraint e = 0.1 as function of the number of RAPs, 
whose signals are centrally processed. Results in the left sub- 
hgure consider a scenario characterized by a Rayleigh fading 
channel with 7 = 10 dB, while in the right sub-figure a 
distance dependent path loss model is taken into account, as 
described in Sec. IllI-BI For both models, three values of A 7 
are considered. For the right sub-hgure, the path loss exponent 
is fixed to p = 2 , the compensation factor is s = 0.1 and the 
average SNR at unit distance is set to 7 ud = 0 dB. In Fig.|5] the 
solid lines are evaluated analytically as derived in Section lTlI-AI 
and IIII-BI while the dots are obtained through simulations 
using the methodology described above. The notches on the 
right side of each sub-figure show the behavior as 00 . 

Fig. |5] shows a good agreement for both models between the 
simulated and the analytical results, emphasizing the validity 
of the proposed analytical framework. 

In Section III-Bl the complexity is assumed to be deter¬ 
ministic for a fixed SNR, while in an actual implementation 
this would be a random variable with expectation C( 7 , A 7 ). 
This additional source of randomness is not considered in 
the analytical framework, and it is one of the reason for the 
slight mismatch between the empirical and analytical results. 
Another reason for the slight mismatch is that a lower-bound 
on the complexity is not applied; an actual decoder always 
performs at least one iteration which determines the minimum 
required complexity. If C( 7 , A 7 ) is a random variable, the 
computational variance would increase. This would conse¬ 
quently cause an increment of the computational gain for a 
centralized RAN. Furthermore, if the lower bound is applied, 
it causes a slight increment of the expected complexity, even 
if the variance might decrease, which lead also in this case to 
a higher computational gain. 

V. Results 

In this section, the analytical framework developed in Sec¬ 
tion m is applied to perform an evaluation of a centralized 
system using some of the metrics introduced in Section III-CI 
In the following, the path-loss exponent is hxed to 7 = 2, the 
compensation factor is s = 0 . 1 , the per-cell constraint on the 
computational outage is set to e = 0.1, and the average SNR 
at unit distance is 7 ud = 0 dB. 

A. Computational Gain 

At first, the computational gain ^comp and its dependence on 
the number of centralized BSs and on the SNR margin A 7 is 
considered. The first sub-figure shows the computational gain 
as function of the number of centralized BSs. The second sub- 
figure shows the computational gain as function of the target 
computational outage probability. 



Nc 

(a) Gain as function of Nc 



logio Ccomp 

(b) Gain as function of e 


Fig. 6. Complexity gain for 7 ud = OdB, rj = 2, s = 0.1. The notches on 
the right side of the first sub-figure show the behavior as TVc —>• oo. 

Fig. | 6 la) shows that the computational gain increases 
sharply for small values of N^, before it becomes flat for larger 
values of N^. Furthermore, it shows a high computational gain 
compared to a distributed system, emphasizing the benefit of 
a centralized solution. The notches on the right side of the 
sub-figure show the behavior as oo. 

Fig. | 6 jb) shows that the computational gain increases as the 
target computational outage probability decreases, since more 
computational resources are required to make sure that also 
peak-requirements are satisfied. This is a consequence of the 
fact that the complexity scales with — logiQ(e) in Q. In this 
case, centralization provides more benefits as diversity effects 
can be exploited to load balance between the individual BSs. 
Furthermore, the computational gain differs more significantly 
for different SNR offsets as the target computational outage 
probability decreases. This effect is mainly due to the fact 
that at lower target outage probability it is more challenging 
to operate close to capacity because the decoder would be 
required to have a steeper error rate curve. 

B. Computational Diversity 

Fig. |7] shows the computational diversity as function of 
the target computational outage probability. For each SNR 
margin A 7 and target computational outage probability, the 
computational diversity is computed by numerically evalu¬ 
ating the derivative in (fTSl i. once the computational outage 
probability is obtained from (l28T l. The computational diversity. 
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logio ^comp 


Fig. 7. Computational diversity as function of the target computational outage 
probability. 
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Nc 

Fig. 8. Complexity-rate tradeoff (CRT) as a function of Nq. The notches on 
the right side of the figure show the behavior as Nc ^ oo. 

similarly to the computational gain, decreases by increasing 
the target computational outage probability. This is due to the 
over-provisioning of resources at lower target computational 
outage probability. Furthermore, the computational diversity 
differs significantly for different SNR offsets as the target 
computational outage probability decreases. 

C. Complexity-Rate Tradeoff 

Fig. H] shows the CRT as function of the number of central¬ 
ized BSs, which quantifies the increase of the achievable rate 
as computational resources are added. Fig. | 8 ] shows also the 
dependence of the CRT on the SNR margin offset. The CRT 
defined in (fl4l i is the ratio between the slope of the average 
rate and the slope of the outage complexity both computed at 
a given SNR margin. In Fig. [ 8 ] the two slopes are individually 
evaluated numerically using (fOT l and (l28T l for each value of 
Nc and SNR margin A 7 . The notches on the right side of the 
figure show the behavior as Nc ^ 00 . 

Fig. | 8 ] shows that the CRT increases as the margin offset is 
increased and this is due to the fact that the decoder operates 
farther away from capacity and therefore does not cause high 
peak complexity. Hence, additional complexity is beneficial in 
order to gain additional achievable rate. In other words, the 
closer the decoder operates to capacity, the more resources 
must be invested to increase the achievable rate. Fig. [ 8 ] points 
out that an optimum point to operate a system exists. Assume 
that a minimum average achievable rate is required. In this 


case, the operator may use a less complex decoder, which 
reduces the required computational resources (and therefore 
reduces costs for information technology (IT) infrastructures). 
However, a minimum computational complexity should be 
maintained since the penalty in terms of achievable rate may 
otherwise be significant. In addition, due to the lower per- 
cell achievable rate, more BSs must be deployed, which 
implies higher deployment costs. Depending on the costs for 
IT infrastructure and mobile network infrastructure, an optimal 
operating point can be determined where deployment costs are 
minimized. 

VI. Conclusions 

Computational limitations in mobile wireless networks have 
a significant impact on the network performance. Such com¬ 
putational limitations are of particular interest to small-cell 
networks, where RAPs are computationally limited due to 
economic constraints. The computational constraints are also 
of fundamental importance to the design of centralized-RAN 
architectures. 

In this paper a new analytical framework has been developed 
to evaluate the data processing requirements for a centralized- 
RAN. Several performance metrics have been introduced to 
help quantifying the benefit of such a network compared 
to a conventional RAN: computational outage probability, 
computational gain, computational diversity and complexity- 
rate tradeoff. The novel concept of computational outage helps 
to quantify the complexity-throughput tradeoff in centralized 
RAN platforms. Such a quantification allows for specifying 
computational requirements for centralized RAN platforms in 
a meaningful way. Furthermore, insights into the computa¬ 
tional gain and computational diversity offered by centralized 
processing are provided. Finally, the complexity-rate tradeoff 
is used to measure the additional complexity needed to further 
improve the achievable rate. 

The analytical framework has been verified by comparing it 
numerically with a large-scale system level simulator compli¬ 
ant with the 3GPP LTE standard. It has been shown that the 
analytical framework provides a good match with simulations 
despite simplifying assumptions which are necessary for the 
tractability of the derivations. 

The analysis in this paper does not account explicitly 
for channel outage. For instance, the required computational 
resources may be higher in the case of a channel outage, since 
in this case the maximum allowed number of iterations needs 
to be used. In a more realistic scenario, the computational 
outage and channel outage can both occur, but only the system 
outage is relevant for the operator. We leave it as future work 
to include channel outage into our analytical framework. Fur¬ 
thermore the effect of co-channel interference could be taken 
into account, e. g., by using tools from stochastic geometry, 
e.g. m and bounds or asymptotical expressions could be 
derived for the novel performance metrics introduced here. 

Appendix A 

This section provides details leading to (fThl l. By knowing 
the PDF of the instantaneous SNR 7 provided by O and 





































































11 


by using the complexity model in Section III-BI the expected 
decoding complexity conditioned on the average SNR 7 can 
be evaluated as follows 


After few algebraic manipulations and using the change of 

ttfcf “t” I—I 

variable 7 =-^—, (1401) becomes 

akl 


00 

— [ ^j( 7 ) ^7)/7i(7)t^7 l3{xi,X2) = - —-77-exp [ -exp{—t)dt— 

7 log ( 2 ) Vafc7/ J _ t 

Nr ^'=+1 oc 

= X! y ^*(7>^7)/7.(7)rf7 (31) f 


-t(xi) 


exp{—t)dt 


Ik 


By substituting (flSl l and (|3l) into (OTT l. the integral in (ISTT l 
becomes 


-r(x2) 
a/cX + bk 


7?+! 


where t{x) = 

The integrals in (14Tb are exponential integrals, defined as 
(i.e, ED) 


J C,(7,A7)/^.(7)d7 = [Ji (7f^^,7«) + 

I 2 (7^11 7 ^)] ■ ( 32 ) By substituting (l 42 b into (l 4 Tb . and (l 4 Tb into (l 39 b yields 


7.« 


El (a;) = 


■ dt. 


(42) 


where 


XI 

h (a:i,X 2 ) = log 2 J y exp dj. 

X2 


Xi 

h{xi,X 2 ) = loga [/fe( 7 )]exp d 7 . (34) 


with 


4(7) = l 0 g 2 (1 -b 7) - {ik^ A7) . 


where 


Ofe = 


4 ( 7 ) ~ Ckin) = akj + bk 


4(7fc+i)-4(7fc) 

_ 07-R 

//c+l 'k 


I 2 (a;i,a; 2 ) ~ 21 og 2 [cfc( 7 )] exp -- 


where 


h{xi,X2) = 2 y |^log2[4(7)]|exp|^-^^ d7. 


4(a;i,X2) ~ 21og2[cfe(7)]exp -- 

V 7 


( 43 ) 


( 35 ) 


Eq. (l 34 b cannot be solved in closed form. However, a piece- 
wise linearization of (iTSb can be done, which leads to 


Finally, an approximated closed form expression for the ex¬ 
pected decoding complexity conditioned over the average SNR 
is obtained by substituting (l32b . (l33b . and (l43b into (OTb . 

Appendix B 

This section provides details leading to (fTTb . The variance 
of the decoding complexity conditioned over the average SNR 
7 is 

Var^^ ,A 7 | 7 {Ez} ®^ 7 ,, A 7 jy *|^ [Ei ^ 7 ,, A 717 {Ei} j ^ 

= ®'^i,A7l7 {Cf } — E^i,A7|7 {Ci} (44) 
where E..^. {Ci} is given by ([16]) and 


( 36 ) 


( 37 ) 


E, 


'li 


bk = Ikil^)- [4(7f+i)-4(7f)]-(38) 

7fe+i Ik 

Substituting (l 36 b into (l 34 b and by integrating by parts yields 


00 

.A7l7K^} = y 

0 

Nr 

= J ( 45 ) 


7fc 


By substituting ([T 5 b and (O into ( l 45 b yields 


-h{xi,X 2 ) ( 39 ) 


E _ 1 (^ 21 '^ 2r(7«,A7) 

k^l 


l0g2 (C - 1) 


14 : {'Jk+n Ik) ( 46 ) 


where 


XI 

h{xi,x2) ==y yos 2 


log 2 (1 -b 7 ) - r ( 7 , A 7 ) 


1^2 


G 


log (2) 7 afe 7 -b 6 fc V 7 / 


exp (—— <77 

7/ 


(47) 


and G = J 
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The integral in ( |47] ) cannot be solved in closed form. 
However, a piece-wise linearization of dTSl l can be done, and 
63 becomes 

Xi 

l 4 (xi,X 2 ) ~ = y (“7) 


By integrating 63 by parts yields 
I4(X1,X2) ^ 

h {xi,X2) 

where 


, 2 f ak bk\ / 7 

-log, 


h{xuX2) = Jexp ^logl (^7+ dj. 


[log(2r 


exp 


7 


ak 


1J akj + bk 


log ( ^7 + ^ ) d'y. 


h{xi) 


h{x2) 

By integrating again by parts, (EB becomes 


I 5 ixi,X2) = 


2 exp 


) 




[log (2)]' 

I<i[h {xi),h{x2)] 

where Ei(-) is an exponential integral given as 63 and 


-dx = -|[ 7 e + log(a 6 )]VC( 2 )} 

i-abr 


+ 


E 

n—l 


7,177,2 


00 . 

71.'n.^ 


n\n 


1 1 1 
2 2 2 


\x 


pFg{ai,... ,ap-,bi,... ,bq; z) is the generalized hypergeomet¬ 
ric function defined as (i.e, ED) 


pXq 


«! 

bi 


:x 


= E 

n =0 


(^l)n • ■ • ^ 

(^1 )n ■ • ■ (^17)71 


(56) 


By applying (l54l l and (ISST i into 63 yields 


76(a;i,a;2) = “2 


7 e + log ^ 

0*7/. 


■C(2) 


(49) 


— 33^3 

akl 


1 1 1 _ 
2 2 2 ’ 


t 

akl 


(57) 


While the generalized hypergeometric function is itself an 
infinite series as well as the Riemann zeta function, they are 
widely known and they are implemented as a single function 
call in most mathematical programming languages, including 
Matlab. 

Finally, by substituting 63 and ( 153 into ( |49] ). and after 
few manipulations yields 


(50) 


I4 {xi,X2) 


[log (2)]^ 


-^ 3 x 3 

akl 


After making the substitution t = akl + bk into (fSOl l and by 
defining the function h (x) = akX + bk yields 

r . ^ 2 fbk 

h{xi,X2) = -5-exp —^ 

[log (2)]2 ^\akl 


(-7 

\akl 




111 ^ t 
222 ’ akl 

t 


log — 

akl 


C(2) 




- [log, {t)f exp 


bk-t 

akl 


h(xi) 

h(x2) 


(58) 


h(xi) 

+ 

h(x2) 

(52) 


Appendix C 

This section provides details leading to ( |23 . The CDF of 
the instantaneous SNR 7 averaged over the spatial distribution 
can be obtained by 


Fj(i) = J F^ (7|a;) fr, (w) dw 


(59) 


where Fj ( 7 |w) = 1 — exp 7 > 0 and 7 is given by 

(I 25 I 1 . Substituting (l26l l into (l59l l yields 


/6(xi,a;,) = (53) 

X2 

The integral (153 can be solved using the equality ll25l Fq. 
(18), Sec. 4.1], 

Fi (ax) . 1 


F^il) = f 

Jo . 


= 1-2 


1 — exp ^ 

—21 

7 ud 

; / exp [ 


^0 \ 

7 ud 


2uj dw 


(54) 


where 7 e is the Fuler-Mascheroni constant and ([ (•) is the 
Riemann zeta function ll23 . Furthermore, let use the identify 


where rj is the path loss exponent, s is the compensation factor 
for fractional power control, and 7 ud is the average SNR at unit 
distance away. 

After few algebraic manipulations and using the change of 
variable t = — (l60l) becomes 

7 ud 

F-yil) = 1 — a ^ [ exp (—f) df— 

V 7ud / J 7 / 7 iid 


(55) 


7ud/ 

poo 

/ exp (—1) f 

Jo 


7/7iid 


where a = 


(61) 


where, by using the Pochhammer symbol (a)o = 1 
and (a)„ = a{a + l)(a -f 2 )...(a -f n — l),n > 1 , 


77(l-s) 


and the integrals in dMI ) are gamma 
incomplete functions. Finally, by using the gamma incomplete 
functions in dMT) . dZTl) is obtained. 
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